SPECTRAL THEORY OF DISCRETE PROCESSES 



PALLE E. T. JORGENSEN AND MYUNG-SIN SONG 



Abstract. We offer a spectral analysis for a class of transfer operators. These 
transfer operators arise for a wide range of stochastic processes, ranging from 
random walks on infinite graphs to the processes that govern signals and re- 
cursive wavelet algorithms; even spectral theory for fractal measures. In each 
case, there is an associated class of harmonic functions which we study. And 
in addition, we study three questions in depth: 

In specific applications, and for a specific stochastic process, how do we 
realize the transfer operator T as an operator in a suitable Hilbert space? And 
how to spectral analyze T once the right Hilbert space Ti. has been selected? 
Finally we characterize the stochastic processes that are governed by a single 
transfer operator. 

In our applications, the particular stochastic process will live on an infinite 
path-space which is realized in turn on a state space S. In the case of random 
walk on graphs G, S will be the set of vertices of G. The Hilbert space Ti. on 
which the transfer operator T acts will then be an space on S, or a Hilbert 
space defined from an energy-quadratic form. 

This circle of problems is both interesting and non-trivial as it turns out 
that T may often be an unbounded linear operator in TH.; but even if it is 
bounded, it is a non-normal operator, so its spectral theory is not amenable to 
an analysis with the use of von Neumann's spectral theorem. While we offer a 
number of applications, we believe that our spectral analysis will have intrinsic 
interest for the theory of operators in Hilbert space. 
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1. Introduction 



In this paper, we consider infinite configurations of vectors (/fe)feez in a Hilbert 
space TL. Since our Hilbert spaces Ti. are typically infinite-dimensional, this can be 
quite complicated, and it will be difficult to make sense of finite and infinite linear 
combinations X^fegz Cfe/fe- 

In case the system [ fk) is orthogonal, the problem is easy, but non-orthogonality 
serves as an encoding of statistical correlations, which in turn motivates our study. 
In applications, a particular system of vectors fk may often be analyzed with the 
use of a single unitary operator U in Ti. This happens if there is a fixed vector 
(fi & TC such that fk = U'^ip for all fc g Z. When this is possible, the spectral 
theorem will then apply to this unitary operator. A key idea in our paper is to 
identify a spectral density function and a transfer operator, both computed directly 
from the pair ((p, U). 

We show that the study of linear expressions Ckfk may be done with the aid 
of the spectral function for a pair {ip, U). A spectral function for a unitary operator 
U is really a system of functions (p^,), one for each cychc subspace H{(p). In each 
cyclic subspace, the function p^p is a complete unitary invariant for U restricted to 
Ti.{(p)'. by this we mean that the function encodes all the spectral data coming 
from the vectors fk = U''(p, fc € Z. For backgrovmd literature on spectral function 
and their applications we refer to [Il[l0l[l6, 19, 20, 21J. 

In summary, the spectral representation theorem is the assertion that commuting 
unitary operators in Hilbert space may be represented as multiplication operators 
in an L^-Hilbert space. The understanding is that this representation is defined as 
a unitary equivalence, and that the L^-Hilbert space to be used allows arbitrary 
measures, and will be a Hilbert space of vector valued functions, see e.g., [6]. 
Because of applications, our systems of vectors will be indexed by an arbitrary 
discrete set rather than merely integers Z. 

We will attack this problem via an isometric embedding of H into and L^-space 
built on infinite parths in such a way that the vectors fk in Ti. transform into a 
system of random variables Zk- Specifically, via certain encodings we build a path- 
space n for the particular problem at hand as well as a path space measure P 
defined on a cr-algebra of subsets of fl. 

If Ti consists of a space of functions / on a state space S, we will need the 
covariance numbers 
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where Z„ : — s- S*, i.e., where the stochastic process is 5*- valued. The set S is 
called the state space. 

The paper is organized as follows. In section [51 for later use, we present our 
path-space approach, and we discuss the path-space measures that we will use in 
computing transitions for stochastic processes. We prove two theorems making the 
connection between our path-space measures on the one hand, and the operator 
theory on the other. Several preliminary results are established proving how the 
transfer operator governs the process and its applications. 

The applications we give in sections [3] and 2] are related. In fact, we unify these 
applications with the use of an encoding map which is also studied in detail. It 
is applied to transitions on certain infinite graphs, to dynamics of (non-invertible) 
endomorphisms (measures on solenoids), to digital filters and their use in wavelets 
and signals, and to harmonic analysis on fractals. 

The remaining sections deal primarily with applications to a sample of concrete 
cases. 

2. Stochastic Processes 

A key tool in our analysis is the construction of path-space measures on infinite 
paths, primarily in the case of discrete paths, but the fundamental ideas are the 
same in the continuous case. Both viewpoints are used in [12j . Readers who wish 
to review the ideas behind there constructions (stochastic processes and consistent 
families of measures) are referred to [H [9l [7] and [18] . 

Let (r2,jF, P) be a Borel probablity space, compact Hausdorff space. (Expec- 
tation E(-) = -dP.) 

Let {Zk)k>o be a stochastic process, and 

(2.1) J^n = a-alg.iZklk < n} 

the corresponding filtration. Let An the subspace in L'^(n,¥) generated by 
JF„. Let Pn be the orthogonal projection of L'^{fl,P) onto An', then the conditional 
expectations E(-|JF„) is simply = P„. 

We say that {Zk)k>o has the generalized Markov property if and only if there 
exists a state space S (also a compact Borel space): 

Zk-.n^ S 

such that for all bounded functions / on S, for all n E N>o, E(/|JF„) = E(/|Z„). 

To make precise the operator theoretic tools going into our construction, we must 
first introduce the ambient Hilbert spaces. We are restricting here to processes, 
so the corresponding stochastic integrals will take values in an ambient L^-space of 
random variables: For our analysis, we must therefore specify a fixed probability 
space, with tr-algebra and probability measure. 

We will have occasion to vary this initial probability space, depending on the 
particular transition operator that governs the process. 

In the most familiar case of Brownian motion, or random walk, the probability 
space amounts to a somewhat standard construction of Wiener and Kolmogorov, 
but here with some modification for our problem at hand: The essential axiom 
in Wiener's case is that all finite samples are jointly Gaussian, but we will drop 
this restriction and consider general stochastic processes, and so we will not make 
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restricting assumptions on the sample distributions and on the underlying proba- 
bility space. For more details, and concrete applications, regarding this stochastic 
approach and its applications, see sections and S] below. 

We begin here with a particular case of a process taking values in the set of 
vertices in a fixed infinite graph G: |13j 

2.1. Starting Assumptions and Constructions. 

(a) G = (Go, Gi) a graph, Go = the set of vertices, G^ — the set of edges. 

(b) {S, Bs, ^j) a probability space. 

(c) The transition matrix is the function 

p{x,y) := P({7 e r!|Z„(7) ^ x, Z„+i(7) - y}) 

defined for all {x,y) G G^, and we assume that it is independent of n. 

(d) From (a) and (b), we construct the path space 

n :={"/= (xoxixa • • • )\ix,^iXi) € G^Wi G N}, 

and the path- measure P = P^. The cylinder sets given by the following data: 
For Ei €Bs, EiCS, set 

P(G(Si,--- ,S„)) 

p{xQ,xi)pixi,X2) ■ ■ ■ p{xn~i,x„)dn{xo)dn{xi) ■ ■■dll{Xn) 



lEo J El 

(e) Starting with {n,T,P), HQ C is a subsigma algebra, let E(-|tJ) be the 
conditional expectation, conditioned by Q. 

If (Xi) is a family of random variables, and Q is the cr-algebra generated by 
(X,) we write E{-\{X,)) in place of E(-|^). 

(f) Let {il,!F,P, {Zn}) be as above. We say that (Z„) is Markov if and only if 

E(/ o Z„+i|{Zo, • • • Zn}) = E(/ o Z„+i|Z„) for all n e Nq. 

(g) From (b) and (d) we define the transfer operator T by 



(2.2) (r/)(x)= / p{x,y)f{y)dfiiy) 

Js 

for measurable functions / on 5*. If 1 denote the constant function 1 on S, 
then Tl^l. 

(h) Let {S, Bs, fi) and T be as in (g), see (|2.12p . A measure /xq on S is said to be 
a Perron- Frohenius measure if and only if 

(2.3) / {Tf ){x)diio{x) — / f{x)diio{x), abbreviated /lo o T — /lo- 
Js Js 

(i) Let {n, !F, P) be as above, and let T be the transfer operator. If fj,Q is a Perron- 
Frobenius measure, let P^'^") be the measure on O determined by using as 
the first factor, i.e., 

p{xQ,Xi)p{xi,X2) ■ ■ ■ p{Xn^i,Xn)dllo{xQ)dll{xi) ■ ■ ■ d^l{Xn) 



Eq <J El 

P,„(G(^i,-- - ,En))d^I{xo). 



In many cases, it is possible to choose specific Perron- Frobenius measures /xq, 
i.e., measures /io satisfying 



^o('S') = f and / {T f){x)d^Q(x) = / /(a;)d/io(x). 
JS JS 

(Note the normaUzation!) 

Theorem 2.1. (D. Ruelle) [2] Suppose there is a norm \\ ■ \\ on bounded measurable 
functions f on S such that the \\ ■ \\-completion L{S) is embedded in L°°{S), and 
that there are constants a € (0, 1), AI G K+ such that 

\\Tf\\<a\\f\\+M\\f\\^, 

where \\ ■ ||oo is the essential supremum-norm. Then T has a Perron- Frobenius 
measure. 

Theorem 2.2. Let (S, fi) be a probability space with S carrying a separate a-algebra 
Bs and /i defined on Bs ■ Let be the path space, and supposed the transfer operator 
T has a Perron- Frobenius measure iiq, then 

(2.4) E(^'')((^oZ„)(^oZ„+i)) = (v9,r7A)i2(^„) 

for all ip,ip ^ L^{fi), and all n G Nq. Here E(i^) :— F{uj)dF(uj) for all integrable 
random variables : — > C; E for expectation. 

Proof. 

E('^")((^oZ„)(^oZ„+i)) 

P{xq, Xi)- ■ ■p{Xn~l,Xn)p{Xn,Xn+lW{.Xn)i>{x„+l)d^lQ{xo)d^Ji(Xl) ■ ■ ■ dpi{Xn+l] 



s J s 



s Js 



p{xo,Xi) ■ ■ ■ p{Xn-l, Xn)(p{Xn){Tlp){Xn)dp.o{xo)dfJ.{xi) ■ ■ ■ dfl{Xn) 

r\Tp ■ iTi;)){xo)d^ioixo) 
'(f{x){Tip){x)d^Q{x) by Perron- Frobenius 



□ 



It is not necessary in (|2.4p to restrict attention to functions iy9,-0 in _L^(/io)- The 
important thing is that the integral (p{x){Tip){x)d^Q{x) exists, and this quantity 



may then be used instead on the RHS in 

Let (Z„)rieNo be a stochastic process, and let Tn be the cr-algebra generated 
by {Zk\ < k < n}. Futhermore, let E(-|J>i) be the conditioned expectation 
conditioned by J-n- 

Theorem 2.3. Let iZn)neNo stochastic process with stationary transitions and 
operator T . Then 

(2.5) E(/oZ„+i|.F„) = (r/)oZ„ 

for all bounded measurable functions f on S , and all n G No 



Proof. We may assume that / is a real valued function on S. Let An ■— all bounded 
J>i-measurable functions. Then the assertion in (|2.5p may be restated as: 



(2.6) / (^(/oZ„+i)dP= / (^((T/) o Z„)(iP 

Jn Jn 

for all ip e An- 

If G An, '^(•) ~ ^{xq, xi, ■ ■ ■ Xn)', and then the LHS in (j2.6[) may be written as 

P{xq,Xi) ■ ■ ■p{Xn,Xn+l)^{xQ ) f {Xn+l)dflo{xo)dfl{xi) ■ ■ ■ dfl{Xn+l) 

p{xo, Xi) ■ ■ ■p{Xn-l,Xn)^{xo ^ Xi ' ' • Xji ){T f ){xn)dfio{xo)dfi{xi) ■ ■■d^{Xn) 



S J s 



ISJS JS 

= / • (T/) o Zn dF. 
Jn 

Hence (123| follows. □ 



Corollary 2.4. Let (11, JF, P, (Zn)) be as in the theorem. Then the process (Zn) is 
Markov. 

Proof. We must show that 

E{foZn+l\Tn)=E{foZn+l\ Zn) 

By the theorem, we only need to show that 

E(/oZ„+i|Z„) - (T/)oZ„. 

In checking this we use the transition operator T. As a result we may now assume 
that (f has the form ip — g o Zn for g a measurable function on S. Hence 



ifiif o Zn+i)dF = I [g o Zn){f o Zn+i)dF = {g, Tf)L2(^) 
g{Tf)dy. = I {go Zn){{Tf) o Z„)dP 



p{{Tf) o Zn)d¥ 

which is the desired conclusion. □ 
Definition 2.5. We say that a measurable function / on S is harmonic if Tf — f . 

Definition 2.6. A sequence of random variables (i^„) is said to be a martingale if 
and only if E(F„+i|J^„) = F„ for all n e Nq. 

Corollary 2.7. Let {Zn)nefio stochastic process with stationary transitions 

and operator T . Let f be a measurable function on S . 

Then f is harmonic if and only if (/ o Zn)ne'Ho is a martingale. 

Proof. This follows from (|2.5p combined with Definition 12.61 □ 

Corollary 2.8. Suppose a process {Zn)neNa stationary with a fixed transition 
operator T : L'^{^i) -> L'^{y). Then ^ ^¥ o Z:^'^ for all n e Nq. 

Proof. Let / and g be a pair of functions on 5* as specified above. Then we showed 
that 

9fd^i = (go Zn){f o Zn 



which is the desired conclusion. □ 



2.2. Martingales and Boundaries. Let G — (G°, G^) be an infinite graph with 
a fixed conductance c, and let the corresponding operators be Ac and Tc- 

Let h : G'^ ^ R is a harmonic function, i.e., Ac/i = 0, or equivalently Tch = h. 

As an apphcation of Corollary 12.71 we may then apply a theorem of J. Doob to 
the associated martingale h o Zn, n G Nq. This means that the sequence {h o Zn) 
will then have P- a. e. limit i.e., 

(2.7) lim h o Z„ = v pointwise P a.e. 

n — ►oo 

The limit function u : ^ M will satisfy v{xqXiX2---) = v(xiX2X3 ■ ■ ■), or 
equivalently, 

(2.8) v^voa. 

The existence of the limit in (|2.7p holds if one or the other of the two conditions 
is satisfied: 

(i) heL°°;oT 

(ii) sup„ /j^ \ho Z„pcflP < oo. 

Proposition 2.9. JTj If h : G'^ R is harmonic and if (i) or (ii) hold, then 

(2.9) h{x) ^ i V dF^ for all x G G", 



where P^, = the measure P conditioned with ^0(7) = 2;. The converse implication 
holds as well. 

Proof. Starting with h harmonic, if the Doob-limit v in (|2.7p exists, then it is clear 
that V satisfies (|2.8I) . By Dominated Convergence, (|2.9p will be satisfied. 

Conversely, suppose some measurable w : ^ R satisfies (|2.8p . and the integral 
in (|2.9p exists then 

{Tch){x) = ^p{x,y)h{y) 

= = x,Zi= 2/)E(z;|Zo(-) = y) 

Y,p{x,y)¥.MZi{-)^y) 

y^x 

= X! y)IE(w|^o x, Zi = y) 

yr^x 

= P4<--)) 

showing that h is harmonic. 
2.3. Solenoids. 



by l lTst 



□ 



Example 2.10. Let S' be a compact Hausdorff space, and cr : S* — > S* a finite-to-one 
endomorphism onto S. Let X„{S) be the corresponding solenoid: 

Xa{S) ^ W S, where No = {0} U N = {0, 1, 2, 3, • • • }, 

riGNo 

(2.10) Xcy{S) = {{xk)kmo\(^{xk+i) = Xk}. 
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One advantage of a choice of solenoid over the initial endoniorphism a : S S 
is that a induces an automorphism a : X^{S) — s- X^{S) as follows: 

a{{xoXiX2 ■ ■ •)) — (o'(xo)a;o2;iX2 ■ ■ with inverse a^^{{xoXiX2 ■■■)) — {X1X2X3 ■ 
Let : S* — > [0, 1] be a Borel measurable function, and set 

(2.11) {Twf){x)= J2 W(y)f(y)^ feB{S),xeS. 



y 



Assume 

(2.12) W{y) = l,yx€S. 

cr(y)=x 

For points x G 5, set D{x) := tt{2/|o'(t/) — x}. A measure /x on 5 is said to be 
strongly invariant if 



f{x)dfi{x) 



y 

a(y)=x 



Lemma 2.11. Assume a measure fj, on S is strongly invariant, and let m be a 
function on S. Set Vf{x) m{x)f{ij{x)). Then the adjoint operator 

V* : L^(/i) ^ L\^,) zs iv*f){x) = E My)f{y) 

V J y 

cr(y)=x 

Proof. See [Til. □ 

Set :— Xa{S) and equip it with the cr-algebra T and the topology which is 
generated by the cylinder sets. 

Set Zk-.n^ S, 

(2.13) Zk{xQXiX2 ■■■)■■= Xk, A: e Nq. 
Let _E C S* be a Borel set, and consider 

(2.14) Z^^{E) = {w e Vt\Zk{uj) e E}. 
Then the cr-algebra on is generated by the sets 

(2.15) Z^^{E) as A; and vary. 
Set 

(2.16) Tn cr-algebra < Zk\k < n > 

where <C • ^ refers to the cr-algebra as specified in (|2.14p . 

In 51 = Xo-(5), consider the following random walk: For points x,y G S, a 
transition x y is possible if and only if cr(y) = x; and in this case the transition 
probability is pw{x, y) := W{y). 

Let /i be a probability measure on S. In we introduce the following Kolmogorov 
measure P :— Pw which is determined on cylinder sets as follows 

(2.17) P(C„) ■.^nC{E^,Ei,E2--- ,E,,)) 

(2.18) =( I ■■■( W{xl)W{x2)■■■W{Xr,)d^Ji{xQ)d^l{xl)■■■d^i{Xn) 



More specifically, P is a measure on infinite paths, and 

(2.19) Cn = {t^ = {tjJouJiuj2 ■ ■ ■)\(t{uJi+i) = uJi, Zk(uj) G Sfe, for Q <k <n}. 

Example 2.12. The following is a solenoid which is used in both number theory 
(the study of algebraic irrational numbers) and in ergodic systems. [4]. For this 
family of examples, the solenoids are associated with specific polynomials p G Z[a;]. 

Let 5 := T" where s G N is fixed; and let p{x) — a^x'^ + aix^~^ + • • • + as; ao 7^ 0, 
be a polynomial, p G Z[a;]. Set 



/ 


ao 








... \ 








ao 

















ao 

















ao 










-a2 -ai) 



Consider the shift a on the infinite torus Jlz'^'* ^ (T'^)^, and set 
Xa := {{zn)nez € (r)^|aoz„+i = FZr,}. 

Then it follows that X„(p) is ct- invariant and closed. As a result, Xcr{p) is a compact 
solenoid. 

3. Graphs 

One additional application of these ideas is to infinite graph systems (G, c) where 
G is a graph and c is a positive conductance function. A comprehensive study of this 
class of examples was carried out in the paper [12]. We will adapt the convention 
from that paper: 

G" ; the set of vertices in G; 

G^ : the set of edges in G; 

and c : G^ ^ M+ the conductance function. 

Assumptions. 

(i) Edge symmetry. lix,y £ G" and (x, y) E G^, then we assume that Cx^y — Cy^x- 
Moreover, (x,y) G G^ (y,x) G G^ 

(ii) Finite neighborhoods. For all x G G", the set Nhh{x) = {y G G°|(x,?/) G G^} 
is finite. 

(iii) No self-loops. If a; G G", then x ^ Nbh{x). 
Convention: If a;, y G G", we write a; ~ y iff (x, y) G G^. 

(iv) Connectedness. For all x,y E G*^ there exists {xi}"^Q C Go such that 
{xi,Xi+i) E , i = 0,1, ■ ■ ■ ,n — 1 xq = X and x„ = y. 

(v) Choice of origin. We select an origin o G G" . 

Definition 3.1. • The Laplace operator A = A^. 

(A/)(x) :=^c,,,(/(x)-/(y)). 

• Hilbert spaces: 
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(i) P{G°): functions / : G" -> C such that ||/||2 = J^xeGO \fi^)\^ < °°- 
Set (/i, /2)2 := ExeGO M^fiix). For every a: € G^, set (5, : G" ^ M, 



1 if = X 
if y 7^ a; 

Note that {6x} is an orthonormal basis (ONB) in P(G''). 
(ii) Ti.E- finite energy functions module constants: 

(3.1) 11/111 = ^EE^-.^i/(^)-/(y)i'- 

all a:~y 

Set 

(3.2) (A, f2)E := ^ E E ^-A7^ - hiy))if2ix) - f2{y)). 

• Dipoles. For all x G there is a unique Vj. ^TCe such that 

{VxJ)E^f{x)-f{o), yfeHE. 

In this case, Vx satisfies Avx = ^ So, and we make the choice Vx{o) ~ 0. 
The function : G° — > M is called a dipole. 

Example 3.2. The dyadic tree. 

• A = the alphabet of two letters, bits {0, 1} ~ Z2. 

• G*^: the set of all finite words in ^ : o = = the empty word, x = 
{aia2 ■ ■ ■ an) € G°, Ui G A, a, word of length n; l{x) = n. 

• G^ :— the edges in the dyadic tree. If x = 0, Nhh{x) = {0, 1} two one-letter 
words. If l{x) = n > 0, X = (0102 • • • a„), Nhh{x) = {(oi • • • a„_i), (xO), (x\)}. 
Set X* :— (ai • • • a„_i). 

• Constant conductance. 

This is the restriction c = 1 on G^ . Then 

(A/)(o) = 2/(o)-/(0)-/(l), and 
(A/)(x) = 3/(a;) - f{x*) - f{xQ) - f{xl), 
if a; e G°, and /(x) > 0. 

• Paths in the tree. If a; = (0102 • • ■ On) G G*', there is a unique path ^(x) 
from to x: the path is 

7(x) {(o, ai), (01,(0102)),- ••((ai-- •a„_i),a::)} 

and consists of n edges. 

• Concatenation of words: For x — (01O2 • • • a„), y = (6162 ' ' ' bm) G G". Set 
z = z(a;y) = (oi • • • o„6i • • • 6,„). 

The dipoles (vx) are indexed by a; G G*' \ (o), and Vx{o) — where o is the chosen 
origin. If G = the tree, then o = = the empty word. 

Lemma 3.3. [12 Let x — (0102 • • • a„), at G A, n — l{x); and y = (6162 • • • Wn), 
bi A, m = l{y). Then 

ify^o 

Vx{y) ■= { 2"-™ • (2™ - 1) - 2^ if m < n 
^ 2^ if m > n 
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(ii) v.,GnE, and||«,||| = f(22»-l). 

(ill) {v.,,Vy)E = |(22-i"('(-)Ky)) _ 1) = #(^(^) n7(y)), for allx.,ye G° \ (o). 

Proof. (i) By the uniqueness in Lemma l3.3[ it is enough to prove that the func- 
tion Vx in (i) satisfies {vx,f)E = f{x) — f{o) for aU / g H_e, and therefore 
also 

(3.3) Aw^r = 5x - So] 

and that (ii)-(iii) hold. 

Specifically, we must prove that 

(A«,)(o) = -1, 

{Avx){x) = 1, and 

(Ai;,)(j/) =0, if 2/^ {0,4. 

Each is a computation: 

(Az;^)(o) = 2vxio) - VxiO) - v^il) 

= - (2 • 2"-i - (2" - 1)) 

= 1 

= So{o). 
And if ?/ 7^ o, but m < n, then 
(Au^)(2/) = 3uj:(y) - Vx{y*) - v^iyO) - Vx{yl) 

= 3 • 2"^™ • (2™ - 1) - 2""''"+^ • (2™~i - 1) - 2 • 2""™^^ • (2"+^ - 1) 
= 0. 

Finally, we compute the case y = x as follows: 

{AVor:){x) = SVxix) ~ Vx{x*) ~ Vx{xQ) ~ Vx{x\) 

= 3 • (2" - 1) - 2 • (2"-i - 1) - 2 • (2" - 1) 
=0-3+2+2= 1 
= bx(x) - bo(x). 

We leave the case m — l{y) > n to the reader, 
(ii) A computation using (|3.ip yields 



x\\E 



\2 _ }_ ^-2n-m^2 
m<n 

1.2- ^^-2" 



2 V 1 

= ^(2^"-l) 

proving (ii). 

(iii) Suppose m = l{y) < n = l{x), x,y E \ (o). From (|3.2p . we see that the 
contribution to {vx,Vy)E only includes words z with l{z) < m. 
The desired conclusion 

(z;„^;,)^ = 2-2™#(7(x)n7(y)) 
11 



follows as in (ii). The possibilities may be illustrated in Figure 1 below. 



□ 



4. Specific Transition Operators 

4.1. Transition on Graphs. Let G = (G''\G^) be a graph with conductance 
function c : G^ —> and transition probabilities 

p{x,y):=^^, y{x,y)GG\ 

C[X) 

Note that c{x)p{x,y) = p{x,y)c{y), which makes the corresponding p-random 
walk reversible. 

Lemma 4.1. Assume that #Nbh{x) < oo for all x G G^. Set 

{Tf){x) ■.= Y.pix,y)f{y), 

and let (Z„) be the random walk on G^ with transition probabilities p{x, y) on edges 
{xy) in G, i.e., 

P({7|Z„(7) = X, Z„+i(7) = y}) = p{x, y) for (xy) € G^ 

Let T be the transition operator, and for (p e 1^{G^), set 

if) ■■= '^(^) 

then for pairs of functions fi and f2 on G°, we have 

Eiifl O Zn) ■ (/2 O Z„+i)) = (T"(/i • T/2)) 

with fi and /2 are restricted to make the last sum convergent. 

Proof. Let /i, /2 be a pair of functions (real valued) on G° such that the point- 
wise product /i • (T/2) is in Z^(G°). Then for n € No, we now compute the Z„- 
expectations: For the P-intcgration on path space fi, we have: 

E((/i ° Z„) ■ (/2 o Z„+i)) 

= / (/i0Z„)-(/2 0Z„+i)dP 

Jn 

= 'Ypi^0,Xl)p{xi,X2)---p{Xn,Xn+l)fl{Xn)f2{Xn+l) 
such that Xi-i^Xi 

= ■ ■ ■ 'Ypi^0,Xl)p{xi,X2) ■ ■ ■ p{Xn-l, Xn)fl{Xn){Tf2){Xn) 

Xo Xl Xji 

= ^ T"(/i.r/2)(a;o) 

xoeG" 
= (T"(/i • T/2)) 

□ 
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Theorem 4.2. Let {G,c) be a graph with conductance c : ^ M+. Assume that 
^Nbh{x) < oo for all x € G° , when Nbh{x) {y e G"\y - x}. Set 

p{x,y) := ^^^/j) and (T f){x) ■.^^p{x,y)f[y). 

Set 

l\G'',fi,) = {/ : GO ^ R\x ^ c{x)f{x) e ?i(G°)}, and (/), ^ c{x)f{x). 

xeGO 

Let F'^'^'^ — p('^'^) be the cylinder path-measure on 

:= {{xaXiX2 ■■ ■)\xi € G", Xi^i Xi,i eN} 

where we use /i^ in the first variable Xg, and counting measure on the remaining 
variables. Then 

^^"^Hifl O Z„) • (/2 O Z„+i)) = (A • Tf2)c 

Proof. 

E<^^^\{fl0Zn)-{f2 0Z,,+ i)) 
= ■ ■ ■ X! c{xo)p{xo, Xi)p{xi, X2) ■ ■ ■ p{Xn,Xn+l)fl{Xn)f2{Xn+l) 

XO XI Xn+1 

such that Xi_i^Xi 

= • ■ - c{xo)p{xo, Xi)p{xi, X2) ■ ■ ■ p{Xn-l,Xn)fl{Xn){Tf2)(Xn) 

Xq Xi Xn 

= 5]c(xo)r"(/i.r/2)(xo) 

XQ 

= (r"(/i-T/2)), 

= {fl-Tf2)c 

In the multiple summations X^xo "^xi ' ^.i^ it is just the first ^^^-summation 
that is possibly infinite; in case the vertex-set G° is infinite. Note that the combined 
summations in the beginning of the proof contribute the integration over the set 
of all infinite paths 7 = {xqXiX2 ■ ■ ■) specified by xq ~ xi, xi '--^ X2, 2:2 ~ 2:3, • • • , 
at each step, moving from xi to the next variable, note that Xi+i ranges over the 
finite set Nbh{xi). For more details on this point, see (|4.ip . below. 

In the last step, we used the following formula which is valid on l^{^c)- 

(4.1) (T^),-(^)„ ^el'ipc). 

We prove ((i7T|) : 

{Tf)c^ c{x)Yp{x,y)ip{y) 
xeG° y~x 

yeG° x^y 

= Y ^(y)c(y) 

□ 
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4.2. Transfer Operators. In section[21 we showed that a stochastic process (Z„)„gNo 
on a probabihty space (il, JF, P) induces a transfer operator T. The derivation of T 
is then essentiaUy canonicaL 

Here, the strategy wih be reversed; but now, starting with T, there is a variety 
of choices of associated processes (Z„)„gNo- 

4.2.1. Setting. Let S* be a compact HausdorfF space. Let {S,B)s,fj) be a Borel 
probabihty measure space, and let p : S x S M>o be a continuous function such 
that 



(4.2) I p{x,y)d^{y) = l fi a.e. x. 
Set 

(4.3) {Tf){x) / p{x,y)dfi{y) for ah / e L°^{S). 



Set 

(4.4) n -.^ Hp = {'^ = {xoXiX2- ■ ■)\x.i <E S, s.t. p{xi^i,Xi) > 0}, 

so an infinite path-space with path transitions governed by te function p. 

Let P = (Pp) be the associated cyhnder measure on ilp as defined in section (5] 
For n G No and 7 = {xqXiX2 • • • ) G rip, set 

(4.5) Z„(7):=a;„; i.e., Z„ : fJp ^ 5 
is an S valued random variable for all n G Nq. 



Theorem 4.3. Let p : S x S R>o be as stated in above. Let T be the 

transfer operator {j.S^ . Then the stochastic process (^n)neNo 14-5^ satisfies 
(4.6) 

IE''')((/i o Z„) • (/2 o Z„+i)) = / (r"(A • Tf2)){x)dy.{x) for all A, A G L°-[S). 



s J s 



s J s 



Proof. The details in the computation for (|4.6p follow those in section [21 but the 
reasoning is now reversed. Indeed, 

E(P)((/l°^n)-(/2 0Z„+i)) 

)fl{Xn)f2{Xn+l)dn{xQ)dfi{xi) ■ ■ ■ d^{Xn+l) 
p{xo,Xi) ■ ■ ■p{Xn-l,Xn)fl{Xn){Tf2){Xn)dpio{xQ)dpi{xi) ■ ■ ■ d^{Xn) 

{T-{fi-Tf2)){x)dy.{x). 

□ 

Definition 4.4. Let T be a transition operator satisfying the conditions (|4.2p and 
(|4.3p . and suppose there is a Perron- Frobenius measure /uq on 5, i.e., 

(4.7) fiooT^fiQ. 

We say that T is ergodic if there is only one probability measure fiQ on [S, Bs) 
which solves (|4.7p . 
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If T is ergodic, and /xq is the (unique) Perron- Frobenius measure, then it follows 
from the Pointwise Ergodic Theorem that for all / e L°°{S), the limit 

(4.8) hm r"(/) = ^i^{J)^ 

n — *oo 

pointwise a.e. exits on S, where 1 denotes the constant function 1 on S. 

Corollary 4.5. Letp, T, S, Bs, A*, and (Zn) satisfy the conditions of the theorem. 
Further assume T is ergodic with Perron- Frobenius measure fj,Q. Then 

(4.9) lim E(f)((/i O Zn) ■ (/2 o Zn+l)) = Mfi ■ Tf2) 
n — >oo 

is satisfied for all /i, /2 G L°^{S). 

Proof. To verify (|4.9p . note that E(p)((/i o Z„) • (/2 o Zn+ij) is already computed 
in (j4.6p in the theorem. 

Since /it is a probability measure, the conclusion (|4.9p now follows from (j4.8p . 
i.e., form an application of the Ergodic Theorem. □ 

4.3. Transition on Solenoids. Let {S,fj,) be a measure space, a : S ^ S an 
endomorphism as specified in section [2l Let O := Xcr{S) be the corresponding 
solenoid. Let : 5 ^ [0, 1] be a function satisfying 

(4.10) W{y) = l; 

y,<T{y)=x 

and let P = ¥f^^a,w be the corresponding path measure. 

Lemma 4.6. For the solenoid set Zn : ^ S, Zn{xo,Xi,X2, ■ ■ ■) = Xn, and 

{Tf){x) = J2y a{y)=x'^iy)fiy)' X ^ S. Suppose T has a Perron- Frobenius 
measure fiQ. Then {Zn)neNo stationary with transition operator T. 

Proof. Let /i, /2 be a pair of functions on S satisfying the conditions listed above. 
For the P-integration on path space n{= X„{S)) we then have: 

E^'^°H(/i o Z„) • (/2 o Z„+i)) 

= / E E E W{xi)W{x2)---WiXn+l)fliXn)f2{ )dfio{xo) 

^ Xl X2 2:71 + 1 

(T(xi)=xo(T(a;2)=a:i (T(x„ + i)=a:„ 
= / E E E W{xi)W(x2)---W(Xn)fl{x,,){Tf2){Xn)dllM 

(t{xi ) — Xq(t{x2)—X1 (T{XrL) — XrL-l 

= I {r\fi-Tf2)){xo)dfio{xo) 
= fio{T^fi-Tf2))^fio{fi-Tf2) 

= {.fi,Tf2)mf,o)- 



□ 
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4.4. Encodings. Let G — {G^ , G^) be a graph where we write for the vertices 
and G^ for the edges. Let S" be a set. We say that G yields an encoding of the 
points in S if there are mappings 

(4.11) T° -.G^ ^ S, onto, and 

(4.12) : G° ^ Functions {S ^ S) 
such that for every e = (x, y) € G^ we have 

(4.13) rO(y) = ri(e)rO(x). 
Examples. G — the binary tree, 

(4.14) 5 = No = {0,l,2,...} 

Finite 

(4.15) ={ ^ a;fc2'=|xfc e {0,1}} 

/c=0 

If n G No is given the finite word (xqX\X2 • ■ • ) in (4.6) is computed from the 
Euchdean algorithm for division with 2. 

Points in G*^ are represented by the empty word o, and by all finite words w = 
(xqX\ ■ ■ ■ Xp). Set 

p 

(4.16) t''(w) =^a;fe2'= = nGNo. 

Starting with w = {xoXi---Xp) G G^ , the three neighbors are (?x;0),(wl), and 
w* :— {xqXi ■ ■ ■ Xp-i) truncation, see Figure 2. 
Set 



(4.17) 




= n i-^ n; see (|4.16p : 

= n 1-^ n + 2P+^; and 



Note that in this example, there is an additional pair of mappaings Nq 



(4.18) 



\a°{n) = 2n 
\(T^{n) = 2n+ 1 
corresponding to the encoding mappings: 

(Jo : {xoxi ■ ■ - Xp) 1-^ {OxqXi ■ ■ ■ Xp) 

(4.19) I 

one step longer 

(Ti : (xqXi ■■■Xp)^ {IxqXi - --Xp) 

Remark 4.7. The same construction works mutatis mutandis with A'^'adic scaling 
rather than the dyadic representation of points in No- Moreover, in the represen- 
tation 

p 

(4.20) r^ = ^XfciV^ 

/c=0 

the choices for Xk may be from any complete set of residues modulo iV, i.e., points 
in No/iV-No, or Z/NZ = the cyclic group of order N. The residues {0, 1, • • • , 7V-1} 
is only one choice of many. 
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Encoding of Z. The representation used in (j4.16p above works for Z as well, but 
with the following modification: 



(4.21) T°(xoa;ia;2 ■ ■ ■ Xp) := -2^ + ^ Xk2^ 

Explanation: 



fc=0 



r°(lll - -1) =-2P + Xk2'' with Xfc = 1, < fc < p 

p+1 times k—0 

= _2P + 2f+i - 1 
= 2f - 1. 

Hence, with this convention we arrive at an encoding of Z. 

Graphs vs Compactification: In the examples, we represent points in the vertex 
sets on a graph G by finite words in a specific finite alphabets. A choice of 
compactification fl of G'^ is the set of infinite paths 7, i.e., 7 — {xoXiX2 • • • ) where 
Xi e G°, and {x^„i,Xi) G for all i e N. 

In each of the examples we present, we build measure P on the compactifica- 
tions Q with use of Kolmogorov's extension principle. This is a projective limit 
construction which proceeds in three steps [11| : 

(i) First specify P only on finite words, i.e., on cylinder sets over G" 

(ii) Check that the prescription of P on cylinders is consistent. 

(iii) With Kolmogorov's theorem than extend P to the Borel cr-algebra of subsets 
in generated by the cylinder-sets [TSl [TT] . 

Definition 4.8. In later applications, the following two cases for P will play a role: 
Consider the subset f^Fin in ^ consisting of paths 7 = {x()XiX2 ■ • • ) which terminate 
in infinite repetitions, i.e., 7 € fipin <^ 3 n such that Xi = Xn ^ i > n. The measure 
P is said to be tight if and only if P(f2Fin) = 1- Alternatively, P(f2Fin) < 1- 

Examples Resumed: Wavelets. We adopt the standard terminology for dyadic wave- 
lets in L^(M), specifically ip for a choice of scaling function; see [11]. Let {ak)kei- 
represent a wavelet filter, i.e., satisfying the following three conditions: 



(4.22) '^akak+21 = ^So,i, 

(4.23) ^afe = l, and 

(4.24) ip{x)=2'^akv{2x-k). 

fcez 

The function is in L^(R) and 

(4.25) / Lp{x)dx = 1 



is a chosen normalization. 

Let ^ be the M— Fourier transform. 

The following result is from [11]. Let := the set of all infinite words, and view 
ri as a compactification of the vertex set G" of all finite dyadic words. 
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Lemma 4.9. For every t G M, there is a measure Pt on Q such that 
(4.26) PtixQXi ■ ■ - Xp) = \(f{t + T°(xoXi • • • Xp))\^ 

where : ^ 1, is the encoding of 14-21^ . 
Lemma 4.10. (See ^.) 

(a) Consider the process (Z„) in (il,Pt) from \4-Zb^ with 

Zn {xoXiX2 ■ ■ ■ ) := Xn G {0, 1}. 

infinite word 

Then there is a transfer operator T such that the process is T -stationary. 

(b) Let 



(4.27) W{e'^) := W{t) 



En Jkt 



where functions W onT are identified with 2TT-periodic functions W on M, and 
where (ofe) is some wavelet filter as in ^J^^^-^JTS^. The transfer operator T 
is then given by 

{Twfm - W{^)f{^-) + + n)f{^- + tt). 

We say that W has scaling- degree 2. 

Following I14-18\ I, let a transition from n to n -\- I be given by a choice of 

xe {0,1}- 

Then 

(4.28) Et{ZnZn+i) ^W{t + xn). 

Proposition 4.11. Let (p E L^(R) satisfying ^72^, and suppose \\(p\\2 < 1- Let 
Q, be the compactification derived from the encoding t'^ of Z in \4-21^ and let t G 
(— tTjTt]. Let Pt be the measure on from lli4-26^ . 

Part L Then the following affirmations are equivalent: 

(a) The translates {</?(• — k)\k G Z} form an orthonormal family in L'^iM). 

(b) The measures Pt are tight measures on for all t. 
M Y.n<,z\v{t + n)\^ ^Ifor alltE^. 

Part IL If the measures Pt are not tight, then the translates {Lp{- — k)}kez still form 
a Parseval frame for the closed subspace V{ip) they span, i.e., we have the identity 

2 



E 



(p(x — k)f{x)dx 



\fix)\''dx for all f eV{ip). 



Proof See HI]- □ 

Definition 4.12. Functions W on (— tt, tt] arising as in (|4.27p for a system of 
wavelet coefficients {ak)kez (|4.24p . are called wavelet filters. A wavlet filter W is 
said to be low-pass if fia :~ 5a, i.e., the Dirac measure at = 0, is a Perron- 
Frobenius measure for Tw- 
in general, if is a Lipschitz function, it is known that T\y has a Perron- 
Frobenius measure [3j. 

18 



Example 4.13. Set 

(4.29) Wf (z) -.= -11 + forz = e'^. 

6 

Then Wp is a wavelet-filter under scaling by 3, but it is not a low-pass filter. 
Indeed, the following scaling law holds for Wp- 

^ Wpiw) = 1, \fz = e'^ e 

w'^—z 

We say that Wf has scaling degree 3. 

It is proved in [5] that Wf induces a wavelet representation on an L^-space 
built from the middlc-third-Canter construction, "Cantor-dust" CD^ in M with 
Hausdorff measure 7i", a = i.e., on i^(Cantor dust, 7i"). 

Cantor Dust CD^ . The points x G CD^ C M. are encoded by 

oo 

X = a_fe3'^ + a_fc+i3'^-i + ' ' ' + «o + ^ 

where A: varies in No, and where aj e {0, 1,2} for j £ 'Z, such that —k < j; but 
where aj attains the value 1 only for at most a finite number of places. 

The Perron-Frobenius measure for is singular with support (/io) = T. 

5. Reprocity Rule for the Spectrum 

In the previous section we saw that a wide class of processes are governed by a 
transfer operator T. If the process in question takes places on a graph G = (G°, G^) 
with conductance c, then harmonic analysis on G is phrased in terms of a Laplace 
operator Ac as follows: 



iAJ){x) = J2 ^(^' - for X e GO. 

Lemma 5.1. Let (G^c) and Ac be as above. Set p(x,y) = £i£i^ for {x,y) G 
and let 



then 

And conversely, 



{AJ){x) = cix){fix)-{TJ)ix)}. 
{TJ){x)^f{x)-^{AJ){x). 



Proof. Left to the reader. □ 

Because of reference to harmonic analysis, we present the results in this section 
in terms of A^, but the lemma makes a translation between Ac and Tc immediate: 
For example, a function / on G° satisfies Ac/ = if and only if Tcf = /. Solution 
/ to either one of these equations are called harmonic. 

Definition 5.2. Let 7i be a Hilbert space, and T) a dense linear subspace. An 
operator A defined on 2? is said to be formally selfadjoint if and only if 

(Aw, v) = (u, Av) 

holds for all u,v E V. 
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A further advantage of Ac over Tc is that Ac is formally selftadjoint, (while Tc 
is not!). 

When we say that Ac is formally selfadjoint, this applies to either one of the two 
Hilbert spaces P{G^), and He the energy Hilbert space. 

In the case of TCe, we take for V the linear span of the family {vx\x G G°} C He', 
see Lemma 15.31 and 15.41 

We continue the setup from the previous section: G = (G°,G^) a fixed graph 
with vertices and edges G^. Let c : ^ M+ be a fixed conductance function. 
Let A = Ac be the Laplace operator. Fix an origin o in G*^, and let {w2:}a;eGo\(o) 
be the system of dipoles. 

Lemma 5.3. 12 (Reproducing Kernel) The system {vx}x£G°\{o) forms a repro- 
ducing kernel in the sense: 

(5.1) {vxJ)e = fix) - f{o) for all f e He, 
where He is the energy Hilbert space. 

Proof. The existence of {fa:} is established with an application of Riesz's lemma: 
If a; € G°, there is a path 7(2:) = — > xi — > • • • ^ x„, e; = (xi x^+i) G G^, 
(generally not unique) such that = and Xn = x. 
By Cauchy-Schwarz, we get 

(5.2) I/(^)-/(«)P<Ei77tII/I" 



e- 



Riesz's lemma applied to He, then yields Bv^ & He such that (|5.ip is satisfied. 
We claim that satisfies the dipole equation 

(5.3) Avx^Sx-60, xeG"\{o). 

This imphes (|5.3p . and if Ah — 0, then :— + h solves (|5.3p as well; and 
vice versa. □ 

Lemma 5.4. '12] LetT^Q := span£{Sx}x£G" , andT^E span^{vx}x^G'>\{o)- 

By "span" we mean finite complex linear combinations, so we consider all finite 
summations 

(5.4) Vo^i^axSx}, and VE^i^bxVx}, 

X X 

where {ox} and {bx} denote finite systems of scalars, Ox, bx G C 

Then A yields a density defined hermitian (i.e., formally selfadjoint) operator in 

each of the Hilbert spaces P{G'^) and He- 
Specifically, T>q is dense in P{G^) and 

(5.5) (w, Aw);2 = (Au, w);2 , VMjWGPo- 

Moreover, V is dense uvHe, and 

(5.6) {u, Av)e = {Au,v)e, VujVG'De 

Proof. The symmetry property (|5.5p is immediate from the definition of A. 

We now prove (|5.6p : Since both sides in (|5.5|) are sesquilinear, it is enough, by 
(15. 4p . to prove 

(5.7) {vx,Avy)E = {Avx,Vy)E, yx,yeG°\io). 
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We have 



{Va,,l^Vy)E = {Vx,6y - 60) E 
by \b.6\ 



= {5y - 5o){x) ~ {5y - <5o)(o) 

by iS.li 

= 4(2/) + l 

= {5x-5o,Vy)E 

by symmetry 

= {l\Vx,Vy)E 

by ]5.3l 

which is the desired eq. ()5.7|) . □ 



5.1. Two Hilbert Spaces. Let G = (G",G) be as above; and let c : R+ 
be a fixed conductance function. Let A and T be the corresponding operators, 
A — Ac the Laplace operator, and 

(5.8) (T/)(x) = /(x)-^(A/)(x), xeG°. 

Pick a fixed o G G*^, and let (w2;)2:gG°\(o) be the corresponding reproducing kernet. 

It is important to understand the two operators in the two Hilbert spaces P(G°) 
and He- By (|5.8p . it is enough to consider just A. 

As an operator in P{G°), the operator A has as its domain 

Vq := all finite linear combinations of {Sx}x£G° 
= span {SxjxeGo; 
while the domain in He is 

Ve ■■= span {vx\x e G" \ (o)}. 

Theorem 5.5. (a) The domains in P and in He: 

(i) 'Dq is a dense subspace in P(G^); and 
(a) De is a dense subspace inHs- 

(Hi) If \j,Nbh{x) < 00 for all x E G^ , then A maps Dq into itself: and Ae 
maps Ve into itself. 

(b) For all vectors T^o, we have: 

(i) 

{f,Aif)i2 = c(a;)|^(a;)p - ^^c(x,y)^(a;)93(?/); 

xeG" x y 

x~y 

(ii) {ip,Aip)i2 > 0; and 
(Hi) {ip,Aip)i2 = {Aif,ip)i2. 

(c) For all vectors ip,^ G Ve, we have: 

ft) 

2 



{^,A^)h^= J2 I(A^)( 
xeG°\{o) 

(ii) {lp,Alp)ue ^ 0; and 
( Hi) {if, AiP)he = ( A(p, ■4^)nE ■ 
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x)\^ 



xeG°\(o) 



Proof. The proof of (b)(ii) is a sequence of steps with repeated apphcation of 
Cauchy-Schwarz's inequahty. The proof of (a)(i) is an apphcation of the last equa- 
tion in the proof of Lemma [531 D 

Remark 5.6. The operator A;2 in P, or A^; in He, may be bounded or un- 
bounded. In all cases Ap is essentially selfadjoint in l"^ [12]; but may have 
defect-subspaces. 

5.2. Dichotomy. 

Remark 5.7. [12] For the graph system (G, c) = {tree, 1) the Laplace opera- 
tor (A,I?o) is bounded and selfadjoint in 1^{G^). For the energy Hilbert space 
Ti.E{tree), (A,Ve) is an unbounded Hcrmitian operator. In fact, A is not essen- 
tially selfadjoint on A; i.e., {A,T>e) has a infinite family of distinct selfadjoint 
extensions in the Hilbert space Ti.E- 

Lemma 5.8. Let TL{-, ■) he a complex Hilbert space, and let V be a dense linear 
subspace in Ti.. 

Let L be a closed Hermitian operator defined on V, i.e., L is linear and satisfies 

(5.9) {u,Lv) ^ {Lu,v) yu,v£V. 
Then the spectrum of A is the closure of the set 

(5.10) NS{L):=I^^^ ueV\io) 

Proof. The Hermitian property (j5.9p implies that the spectrum of L is contained 
in M. 

Now suppose Ao G M, and that 

(5.11) dist(Ao,7V5(i)) = ei > 0. 
We will show that A must then be in 

M \ spec{L) — the complement of the spectrum 
= the resolvent set. 

Let u € X'\ (o). Then 

||Aoii - Luf = \l\\uf - 2Xo{u, Lu) + \\Luf 
Setting xi := -^jf^^ G NS{L), we get 

(5.12) IIAo^. - Lu\\' = \\uf ■ (Ao - xi)' - \\u\\^xl + \\Lu\\' 

(5.13) > 2^1,^^1,2 _(!^ 

(5.14) >\\ur-ei 

where we used Schwarz' inequality in the last step; viz., 

or 

Ikll 

By virtue of the inequality (2.11), we may define an operator 
i?o = -R(Ao) : range{\{)I — L) — > H 
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by 

(5.15) Rq{Xu-Lu)^u. 
Extend i?o by setting it = on the ortho-complement 

(5.16) {range{XoI - L))-^ ^ N{Xo- L*). 

Here L* denotes the adjoint operator. 

From (|5.15p . we calculate that Rq : H H defines a bounded inverse to XqI— L, 
and so Aq € resolvent(L); and conversely. □ 

Let {vx}xeG°\(o) be the system of dipoles, and set 

(5.17) M ■.^{{vx,Vy)E) 

viewed as a Hermitian matrix, x — row index, y — column index. 

U^={Cx) (ZP{G°), set 

(5.18) (MO. = 

y 

matrix multiplication, where 

Mx,y := {vx,Vy)E- 
Then M is a density defined Hermitian operator in P{G^). 

Theorem 5.9. Let (G, ^) he given and let A he the corresponding density defined 
Hermitian operator inTiE- Then 

(5.19) spec„^ ( A) C [0, oo) 
and 

(5.20) speoH,{A) ^ {speci.{M)r' 
where we use the charactors — oo, and — = 0. 

^ oo 

Moreover, 

(5.21) (specpiM))-^ = {l/X\X G specpiM)}. 
Proof. For {C^) G J^, set 

(5.22) 

x<£G''\{0} 

Then u £ V, and 

(5.23) {u,Au)nE^YYlx£.y{'"x,Avy)E 

X y 



(5.25) = E 1^- 



2 

> 



Since vectors mTH-E are equivalence classes modulo the constant function on G°, we 
may add the restriction = in (|5.22p . and the operator A will be unchanged. 

The modified equation {2.22) then needs 

(5.26) {u,Au)E^m\l 
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Claim 5.10. 

(5.27) ||u||2^^ = (^,MC)p 

Proof, (of Claim 2.6). We compute: 



X y 

= {(.M()p, 

as claimed. □ 
The desired conclusion (|5.20p now follows: If it G V \ (o) is given by (|5.22p . then 



^'■''^ hill - {i,MO- 

By taking closure, we obtain the sets on the two sides in (I5.20p □ 

Corollary 5.11. If £, — {£,x) G J-{G^ \ (o)), then the representation 

(5.29) u = Y,^xVx 

is unique; in particular, the system {vx)x£G°\{o) linearly independent. 

Proof. Let m S V have a representation (|5.29p as a finite summation with e C. 
Let y€G°\(o). Then 

{Sy,u)E = ^£.x{Sy,Vx)E 

X 

= iv 

In particular, if u = 0, then = 0, Vy G G" \ (o). □ 

Corollary 5.12. If F <Z \ (o) is a finite subset, then is not in the spectrum 
of the matrix 

(5.30) Mp := i{vx,Vy)E)x.,yeF- 

Suppose S spec(Mi?) where F is a fixed as in the statement of the Corollary 
\EM Then 

3^ef{F\(o)) 

such that 

(5.31) (A/0. 

yeF 

Setting u J2yeF ^v'^v note that 

(5.32) u e {{v,}xeF)^ 
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Claim 5.13. 

(5.33) u e i{v^}^<.GO\{o))^- 

We need to prove this only if x e \ i^. 
Combining (|5.27p and (|5.3ip . we get 

\\u\\1 = {^,mOp 

so M = a constant function on G°, and (|5.33p is satisfied. 

6. The Energy-Inner Product 
(G, c) = (treeT, 1), o = 0, c = I. Explicitly form for v^, x e G" \ (o). Set 

(6.1) a; = (010203 ••• a„) e G" \ (o) 0^6/1 = 0,1. 

(6.2) 7(2;) = {(oai), (0102), (0203), • • • , (a„_2an-i), (an-ia™)} 
where 7(0;) is a path. Note 7(x) C G^ = edges in T. 

Example 6.1. x = 101 vertex, {{(p, 1), (1, 10), (10, 101)} = 7(2;) tt7(x) = 3. 
Theorem 6.2. Let (T, 1) be as usual, o = 0, and let He = the energy span 



(6.3) ii/iil-iEE(/(^)-/(y))' 

a: y 

but with edges {x,xb) = e, x E , b E A~ {0, 1}, c(e) = 1. Then the function 

(6.4) v^{y)■.^^{^{x)r^^{y)) 
solves 

(6.5) (i;„/)B = /(a:)-/(o), V/ € T^ij 

(6.6) Av^ = - <5o, X e G" \ (o) 
and 

(6.7) (i;,,«,)i5 = tt(7(a^)n7(2/)) Va;,yeG°\(o) 

Proof. Proof of g^j. By (HH) a; = (0103 . . . a„) G G" \ (o). Let 2: be as in 
Set 7(0;) = RHS in ^ C G^. Neighbors of 

X — > fli . . . a„_i 
— > xO 
— > xl 

If X = a, n = 1, Nbh(x) = {o, aO, al} 
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Cases. 

n=l See Figure 4 

(Aw^)(o) = 2vx{o) - Vx{Q) - Vx{l) 
= 0-1 

= (5j;(o) - 5o{o) 

{ls.Vx){x) = 3wa;(a;) - v^io) - Vx{xQ) - Vx{xl) 
= 3-0-1-1=1 

Msel l6.4t 

= {Sx - (5o)(a;) 

Now, let y £ G° \ {o, x}. y — {bib2 • • • 6fe), 6i G A = {0, 1}. Suppose x <^ y 
Av,{y) = 3-1-1-1 = 



More cases are = 0. 
n>l X = (ai02 • • • a„). A computation yields 

Ax{o) = - vx{0) - v^l) = - 1 = -1 

Ax{x) = 3n - (n - 1) - 2n = 1 
= [Sx - So){x) 

Ax{y) = yeG"\{o,x} 
Several cases e.g. y < x, etc. 

Ax{y) = 3vx{y) - Vxih ■ ■ ■ bk-i) - Vx{yO) - Vx{yl) 
= 3fc - (fc - 1) - (fc + 1) - fc = 0, etc. 
Computation of 

\\Vx\\l = E{Vx) 



s t 



U2;(a;) - Vxio) = n - 
tt(7(^)) 



{Vx,Vy)E = {Vx,AVy)l2 

all finite functions, see ( |d.41 i 

, K'-^a - So) 12 

byi 



= vx{y) - Vx{o) = 

bv i6A[ 

= tt(7(a;)n7(y)) 
b ml6.4[ 



□ 
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Set M = {{v^,Vy)E) = mi{x) n 7(2/))), x,y e \ (o). Given 

specp(M) = (spec„^(A))-i^^^^^^^_^ 
From our theorem above A (unbounded spectrum), closure A = V, V C T-Le- 
Corollary 6.3. Ve C G° \ (o) finite, 3A e specpiMp) such that A < e. 
Note 

Mp = {{Vx,Vy)E)x,yeF 

= (tt(7(x)n7(j/))).,j,eF 

and <^ spcc;2 (Mi?). 

Problem: Find a systematic way of selecting F. See Figure 8. 
It is much easier to find Mp with specpMi? — > 00. 



Example 6.4. 

^1 1> 

(6.8) I 2 I or 



1 1 
1 n 



,10 2 
,+ l± V(«+l)2-4(n-l) 



, and 



_ n + 1 - V(n + 1)2 - 4(n - 1) 
A„ - ^ 

2(n - 1) 



n + 1 - Vn^ - 2n + 2 
1 



Actually both expand part of specpM as intervals. 

7. Karhunen-Loeve 

Definition 7.1. C G° \ (o) F finite, G^ infinite, (G, c) fixed. Fix o G G^ A = 
Ac, He energy Hilbert space (wa,, f)E = f{x)- f{o), gHe- 

Definition 7.2. He{F) := span^^F{''^x} 

General (G, c) — > Fix point o e G° — > A^, He, G° infinite. Fix v^, for a; e 
G° \ (o), determined from Riesz applied to He- 

(7.1) {v.J)E = m-fio), xeG°\{o), 
and consider the infinite matrix 

(7.2) M = {{v,,Vy)E), ar,yeGO\(o) 
and its finite F x F submatrices 

(7.3) MF = {{v^,Vy)E), x,yeF 
so the matrices are 00 x 00, or |F| x 
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Important Formula. Observe 



(7.4) {Vx,Vy)E = Vy{x) - Vy{0) 

(7.5) = Vy{x) = Va:iy); 

in other words kE{x,y) := Vx{y) is a reproducing kernel. 

Since x,y ^ G'^ \ (o); and u^r : — > M (i.e., real valued) convention: Vx{o) = 0. 
Diagonalization motivated by the classical Karhunen-Loeve theorem, see 14) and 

m- 

7.1. Finite-dimensional Approximation. Apply the Spectral Theorem to M 
and Mf- The Hilbert space is P{G") or P{F) ~ C'^l with (^,77)3 = T,x^Vx as 
inner product. 

For {Mp, P{F)) the spectrum is always discrete, and for some cases i.e., (M, P(G")) 
it may not be discrete. 

In the discrete case, there exists M = Mf ONB e ?^(G") or l'^{F) 

eigenvectors 



X 

^ = G ;2(i^) such that 
(7.7) M^^j = AjCj, Ai>A2>--->0 e Z^(-F),||0ll2 = 1 (in the F-case) 

In the infinite case spec(M) for /^(G") may accumulate both at and at 00. 

Since = {vx,Vy)EA € M, we may take all ^fe : G° — > M real valued. Fix 
FcGO\(0): (^l\F). Set 




Lemma 7.3. //i^ is fixed then G P{F) is an ONB. Set 
(7.9) Mf^^ = \U^, 



then 

zs an extension of : F M. from F to G^ . 
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if j = fc 

1 ifj^k 



Proof. By ^ if z e F: 



w^fe (2) = ^^(^-)^kix) 

^ xeF 



by(4.4) A, 



= ^Cf (-) 



1 



□ 



Lemma 7.4. Fia; F C G" \ (0) /inite, and let 

(7.10) i«r(-) = l^I]ef(a:)t'.(-), (1,2,--- ,|F|) as tn Lemma\7;^ 

T/ien {w^jfc is an orthonormal system in He (thus in each of the Hilbert spaces) 
i.e., with the inner product 

(7.11) {u, v)e := ^ XI XI Cxy(w(a;) - u{y)){v{x) - v{y)). 



2 

all xy xr^y 



We have 



Proof. We have: 



1 



byElAj""'" [0 iffc^j 



□ 



Set uj^ — y/XjWj'] then 
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7.2. Normalization. The following different normalization wj^ = y^wj^ satisfies 



(7.13) 

so 

(7.14) 

Note that the 
(7.15) 



lF = V^,e,(-) onF. 



7.3. Projection Valued Measures. Set 

(7.16) (Xj) :^ \uj >< uj\; Dirac notation for rank-one projection, 

so a projection in He on the one-dimensional subspace Cuj^. Then P^{-) is an 
orthonormal projection system, and it has a limit as ^ oo which is a global 
spectral measure. 
We claim that 

(7.17) sa(< ) = {uf, Auf) € specHA^v) 
Lemma 7.5. (Spectral Reprocity) 

1 



(7.18) 
Proof. 



sa('", ) 



1 



xeF 



^ E E A^y)-B 

xGF yeF 



^EE^^-(^)^^(2/)('^-(2') + i) 
E^^(^) 
E^.(- 



xeF yeF 



ml 



b» l7J3l 



Example 7.6. 



= (J 3), sp = {l,3}. 



same spectrum, but different . 



(a> - (6) = 1- 



□ 
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Set Rf{X) := |(^f)|^). Then (ua,Aua) = Rf{X); see Lemma O 

In the examples: 

= I ? I ^^^^^ ^ ] smaller for off-diagonal. 



- ( ^ ) , A± = 5 ± Vl3 



= ( : : I . A 



3 7 



3 I\ , 7±V5 



1 4/' "± 2 



A = 5±yi3^ i?(A) = \ 



A 1 

< 



A l + (2^)2 A 



I TO 

In both cases, we have: 

1 A 



- A ' I + (A-1)2 
A. A_ 



Rf{X-) = — ^ + , /' 7T^ ~ 4 as TO ^ 00. 

TO — 1 1 + (A_ — 1)^ 

We now illustrate by an example that points in the spectrum can go into 00: 

(u,Au) 

If A e speci2{MP) set = ^ J^.^p ^xix)v,{-) M^x = A^a, H^lb = 1 ^ \\ux\\e = 
1 so sa(m) = {u,Au) = i(l + \\Pxe\\l), e = ep = Xf(-), -Pag = {^x,e)2^x 

Theorem 7.7. T/ie truncated operators P-^^(^)Ax)£P-h^(f) ^i^s spectral growth 
~ ©(tJi^); so As is unbounded inTi.E- 

Proof. The idea is to perform a diagonalization of an infinite matrix (M^ y) x,y d 
G° \ (o) ; a method inspired by Karhunen-Loeve [TH |T7] . Here F C G°\ (a) is fixed 
and finite. The following computations refer to F: {£^k) is an ONB in l'^{F) satisfying 
(!7.19p below; set Wk = ^ J2xeF ^k{x)vx, and Vk = V^kWk = -^J^xeF ^kix)v^. 
Then 

(7.19) M^a-A^a, and {Cj,^k)PiF) = S^.k- 

We may now compute the matrices: 

{uj,Auk)E ^^=^^^j(a;)Cfc(x)(w^, Au3,)_E 



Set ExG-F 0(2-') = e)2 = (0) where e = = xf- 
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Then the matrix entries are: Off-Diagonal: 

1 



and Diagonal: 



We further used the following identity: 

FxF 

= S,.k bydZHl). 
This may be summarized in the following matrix form: 



MM + W 



V Ai A3 



/ A2A3 



^ x:(i + (e3)^) 



v 



7 



If for some S E M+, Xj > 6, i.e., bounded from below, then the operator 








• 


■ ^\ 







• 


■ 





1 

A3 


• 


• 


[■ 






■ 0) 



is bounded. So 
(7.20) 



1 



=(e,>(a> 



must be unbounded, i.e., || • ||;2(^)^;2('^-) — )■ 00. But (j7.20p is a rank-one operator; 
|p><p|, p = FcG°\(o) is fixed where p=(ej)e;2(l,2,--- ,tti^), 

i.e., p = and pj" Aj = . 

Now, 



iipn!?.(,...,.) = E%- 



So in conclusion 



lim > —4 — r = CXD 

Pick (5 G K+ and assume > (5. Then we need 



hm E(^f)'--- 

F— >oo 
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We have (j) = 0, M^^f = A/^f , U^h = L (CD = E.eF?r(^)> and 
Ea(^a )' = ttf'; so indeed 

^lim^^(ef)2 = limttf^ = oo. 
Conclusion: spec7^g(jv)(Ai.E) ~ ^F) oo 

□ 
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